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Abstract 



The test of the crab waist collision scheme, undergoing at the e + e~ Frascati DA$NE 
accelerator complex since February 2008, requires a fast and accurate measurement 
of the absolute luminosity, as well as a full characterization of the background 
conditions. Three different monitors, a Bhabha calorimeter, a Bhabha GEM tracker 
and a gamma bremsstrahlung proportional counter have been designed, tested and 
installed around the interaction point end of 2007-beginning of 2008. In this paper, 
we describe these detectors and present their performances in various operation 
conditions during the 2008 and 2009 DA<I>NE runs. 
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34 1 Introduction 



35 Proposals of future flavor factories [Tf2f3] emphasize the need of very high 

36 luminosities. For instance, the new generation of B-factories [T|2] requires im- 

37 provements up to two orders of magnitude above the performances of the 

38 PEP-II [1] and KEK-B [5] e + e~ colliders. Among the ideas currently being 

39 developed to achieve this ambitious physics-driven goal, the crab waist com- 

40 pensation scheme associated with large Piwinski angle and low vertical beta 

41 function [B] is very promising. Luminosities as high as 10 36 cm~ 2 s~ 1 could be 

42 reached with beam currents similar to those operated routinely in today's 

43 accelerators, which would also help keeping the background under control. 

44 In addition to being based on existing technologies, this scheme would sig- 

45 nificantly limit the power (and hence the cost) needed to run such a new 

46 machine. 

47 The DA$NE accelerator, located in the National Laboratory of Frascati 

48 (INFN) and optimized for the production of mesons (-y/s =1020 MeV) at a 

49 high rate, has been modified in 2007 to test the crab waist sextupole compen- 

50 sation scheme [6] with the aim of reaching higher luminosity while controlling 

51 the background. After completion of this upgrade, operations restarted during 

52 winter 2007-2008. 

53 Since 2000, DA$NE has been delivering e + e~ collisions to three experiments 

54 KLOE [7j, FINUDA [8j and DEAR [9], steadily improving performances in 

55 terms of luminosity, beam lifetimes and background. The best peak lumi- 

56 nosity was ~ 1.5 x 10 32 cm _2 s _1 with typical daily integrated luminosities of 

57 ~ 8 pb _1 during the KLOE run. According to calculations [6], the modified in- 
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58 teraction region should increase the luminosity by a factor 3-5 with respect to 

59 the previous colliding scheme. To identify such a significant gain, a luminosity 
eo measurement precise at the ~ 10% level is enough. On the other hand, real 
ei time and accurate information is required regardless of the machine conditions 

62 to allow the DA$NE operators to tune this new collider and to get relevant 

63 feedback for their optimization studies. Some redundancy between the various 

64 measurements is also important, in order to fight transient backgrounds which 

65 could impact strongly a particular detector. 

ee Parallel to the upgrade of the DA$NE interaction point 1 (IP), the SID- 

67 DHARTA experiment (Silicon Drift Detector for Hadronic Atom Research by 

es Timing Application) [10] aiming at studying kaonic hydrogen and kaonic deu- 

69 terium has been installed at the same location. The presence of this additional 

70 detector, whose operation requires a very good shielding against machine back- 

71 ground, has consequences on the design and on the performances of the main 

72 luminometer. These are discussed in the following sections of this article. In 

73 principle, SIDDHARTA can also provide a luminosity measurement based on 

74 the counting of charged Kaon pairs produced by the decay — > K + K~ whose 

75 branching fraction is well-known. However, this method suffers from a few 

76 limitations: a low rate (around 25 Hz at a luminosity of 10 32 cm _2 s~ 1 ); the de- 

77 pendence on the exact center of mass (CM) energy because of the resonance 

78 lineshape; the need to separate efficiently true Kaons from minimum-ionizing 

79 particles, a background strongly machine-dependent; a limited duty cycle for 
so technical reasons. 

si Therefore, various independent luminosity-oriented detectors have been built 

82 around the IP and put in operation beginning of February 2008 with a three- 

83 fold goal: to guarantee an accurate measurement of the absolute luminosity, 
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B4 to monitor the background levels and to provide powerful and fast diagnosis 

85 tools to improve the machine. In the following, we describe the design, the 

ae construction and the operation of these detectors, such as the results achieved 

87 during the 2008 and 2009 DA$NE runs. 

as The paper is organized as follows. Section [2] gives an overview of the processes 

eg used to measure the luminosity and presents the experimental setup. Then, 

90 Sections [3] and H] describe in details the detectors and the trigger /data ac- 

91 quisition systems. The recorded Bhabha events are characterized in Section [5] 

92 which also explains the procedure used to estimate and subtract online the 

93 background. Section O deals with the Monte-Carlo simulation which is crucial 

94 to convert Bhabha event rate into absolute machine luminosity. Validation 

95 studies based on data-Monte-Carlo comparisons are provided in Section [7] 

96 where the systematic error to the luminosity measurements is also computed. 

97 Experimental results are given in Section [8] before we finally conclude in Sec- 

98 t ion [91 

99 2 Overview 

100 2.1 Physical processes used for luminosity and background measurements 

101 Since the early days of e + e~ colliders, well-known electromagnetic processes, 

102 such as Bhabha scattering [11] or single and double e + e~ bremsstrahlung |12j . 

103 have been used to monitor and measure collider luminosity. The Bhabha elastic 

104 scattering e + e~ — > e + e~ has a very clean signature of two back-to-back tracks, 

105 of energy equal or close to the beam energy (510 MeV at DA$NE ). The cross 
loe section o Bhabha of this process has a very steep dependence on the polar angle 
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107 9: UBhabha oc I/O 3 - Most Bhabha events are produced at small polar angle, 

108 an area where focusing quadrupoles or other machine elements are generally 

109 located, which limits the measured event rate. Nevertheless, Bhabha scattering 
no at large angle still produces a sizable counting rate which can be exploited to 
in measure luminosity; considering a polar angle range of 18°-27°, the Born cross- 
n2 section of the process is as high as 5 /zb. 

n3 The e + e~ — > e + e~7 (single bremsstrahlung) process has a cross section of 

H4 169 mb for photon energies above 10 MeV, thus resulting in a very high count- 

H5 ing rate. The photon undergoes a small angular deviation from the charged 

lie track original direction: 95% of the signal is contained in a cone of 1.7 mrad 

H7 aperture. Moreover, its cross section depends only logarithmically on the CM 

us energy and is therefore relatively independent of the actual machine param- 

n9 eters. On the other hand, it suffers heavily from background caused by par- 

120 tides interacting with the residual gas in the beam-pipe or lost by Touschek 

121 effect [T3] . 

122 2.2 Background from Touschek effect 

123 Achieving high luminosity is not enough: future uses of the crab waist scheme 

124 will depend on how well the backgrounds can be controlled in this configura- 

125 tion. For the DA$NE collider, both the experimental machine- induced back- 

126 grounds and the beam lifetimes are dominated by the Touschek effect due to 

127 the use of dense beams at relatively low energy. The Coulomb scattering of 

128 charged particles in a stored bunch induces energy exchange between trans- 

129 verse and longitudinal motions. Small transverse momentum fluctuations lead 

130 to larger longitudinal variations as the effect is amplified by the Lorentz factor. 
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131 Off-momentum particles can then exceed the momentum acceptance given by 

132 the radio- frequency (RF) bucket, or they may hit the aperture when displaced 

133 by dispersion. In both cases they get lost. 

134 Several studies have been performed to control and reduce the background 

135 induced by the Touschek effect and to optimize the signal to noise ratio for all 

136 the experiments that have been running at DA<3>NE. This has been achieved 

137 by adjusting optical parameters like the orbits at the IP or the strength of 

138 the machine sextupoles and by minimizing the radial beam size upstream 

139 from the interaction region (IR). Moreover, the insertion of a proper set of 
wo collimators (see Fig. [Q, together with simulation-based tracking studies of 

141 Touschek scattered particles [H], have also helped reducing the background 

142 significantly. 

143 Particle losses due to Touschek effect at DA$NE are expected to be higher for 

144 the crab waist configuration than for previous optics. However, the collimators 

145 are efficient even if longer jaws would have been more useful around the IR due 

146 to the low emittance. This hardware modification is planned in the near future 

147 for the next KLOE run. In the present setup, a very careful shielding has been 
us designed to allow clean luminosity measurements and to increase the signal 

149 to noise ratio for the SIDDHARTA detector. In addition, lead shieldings have 

150 been added behind the IR collimators to stop as many secondaries as possible. 

151 2.3 Detectors for luminosity and background measurements 

152 In the new DA$NE interaction scheme, the two 510 MeV beams collide with 

153 a total crossing angle of 50 mrad (25 mrad per beam). At the IP, both beams 
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Fig. 1. Layout of the DA&NE main rings with the collimator locations shown by 
red and blue arrows for positrons and electrons, respectively. The IR collimators are 
located upstream of the IP while the other ones are in dispersive regions close to 
wigglers. 

154 are pointing toward the inner part of the ring; they are separated at 57 cm 

155 from the IP. Two low-/? permanent quadrupole magnets, the QDOs, are lo- 

156 cated on each side of the IP before the common beam pipe gets split into 

157 two separate sections. The QDO magnets are approximately cylindrical, with 

158 a length of 24 cm and a 10 cm radius. Their front faces are located at 32.5 cm 

159 from the IP. The SIDDHARTA detector is installed in the space between these 
wo two magnets; its shape is asymmetrical because of physics- and construction- 
al related constraints - see Ref. [10J for details. It mostly covers all the solid 

162 angle above a given polar angle (around 30 degrees) whose actual value de- 

163 pends on its shielding. An overview of the upgraded DA$NE IR including the 

164 luminosity monitors (a Bhabha monitor and two forward gamma detectors) 

165 and the SIDDHARTA detector is shown in Fig. [2j 
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Gamma Monitor 




Fig. 2. Overview of the upgraded DA&NE interaction region, showing the Bhabha 
calorimeter, the GEM tracker and the gamma monitors. The bottom picture was 
taken in July 2008 during an access to the IP; on its top right corner one can see 
the system of coordinates which is used in the following. 



lee Bhabha events are counted by a double-arm lead-scintillator calorimeter in- 

167 stalled around the QDO magnets. This detector can operate in coincidence 

168 with a GEM tracker located closer to the IP to help defining its polar accep- 

169 tance. The luminosity number C meas is extracted in real time by comparing 

170 the background-subtracted Bhabha event rate with the prediction from the 
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171 Monte-Carlo simulation - see Sections [5] and [6] for details. One has 



*~"meas — j^MC ' V ' 



173 where Rmeas is the background-subtracted Bhabha rate measured at the 

174 DA$NE IP and Rgp is the rate predicted by the Monte-Carlo simulation as- 

175 suming a luminosity of C° = 10 32 cm _2 s^ 1 . The reference Monte- Carlo-based 

176 rate includes all acceptance effects: detector geometry, trigger, quality of the 

177 energy reconstruction, etc. 

178 This luminosity measurement can be compared online with a "geometrical" 

179 luminosity estimate, C geom , which is based on beam sizes (as measured by the 

180 synchrotron light monitor, SLM) and beam currents. 

7+ J- 120 1 
isi i~, geom — - — — — - , x ocrabI/J 
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182 with the Ccrab reduction factor given by 



Ccrab = 0.55981 - 0.007474 x 7 med - 8.121 x 10~ 5 x 7 2 ed 
1 ( 1+ r\ 



'med 



2 \ N~ 



183 j3y LM is measured with the quadrupole oscillation method [15]; f3y P is es- 

184 timated using the MAD machine model [16]; the Oy LM factors are com- 

185 puted from the SLM readouts. A^ is the number of bunches in collision (i.e. 

186 min(N + ,N~)) and fnF the RF system frequency, which is divided by 120, 

187 the harmonic number of the accelerator. The coefficients of the Ccrab reduc- 

188 tion factor have been computed numerically; they take into account various 

189 effects (hourglass, non-zero crossing angle, crab sextupoles) and the formula 
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190 is adapted to the DA$NE specific configuration of beams with different beta 

191 functions, emittances and bunch lengths [17]. Data are updated every 15 sec- 

192 onds via the DA$NE control system. 

193 In addition to the Bhabha luminometer, two gamma monitors (PbWCU crys- 

194 tals) are located 170 cm away from the IP on either side and measure the rate 

195 of single bremsstrahlung events. These detectors provide a fast feedback when 
we the machine conditions change; they have also been used to scan the beam 

197 profiles or to estimate the level of Touschek background. A detailed descrip- 

198 tion of these three systems (the Bhabha calorimeter, the GEM trackers and 

199 the gamma monitors) can be found in the following section. 

200 3 The experimental setup 

201 3. 1 Calorimeters 

202 The Bhabha calorimeter is divided into four modules (two on each side of the 

203 IP) which surround the permanent QD0 magnets - see Fig. [2j Its acceptance 

204 in polar angle starts at 18° and goes up to about 27°; the angular upper 

205 bound depends on the exact SIDDHARTA shielding configuration but the 

206 shower containment is degrading quickly above this angle. The modules are 

207 segmented in five azimuthal sectors, each 30° wide. 

208 The calorimeter design is driven by the need of having the main IR platform 

209 supporting its weight (about 400 kg in total) rather than the very fragile 

210 vacuum beam pipe. Each arm is made of two pieces (top and bottom halves) 

211 with a light aluminum skeleton sandwiching a stronger structure located in the 
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212 machine mid-plane. A fan-shaped hole leaves a large and well-defined fraction 

213 of the 180° acceptance free of any passive material. To keep the mechanics 

214 simple, the acceptance is covered by trapezoid (instead of ring-shaped) sectors. 

215 In the chosen design, the supporting structure is made of two parts connected 

216 by a longitudinal bar, one of each side of the fan-shaped hole. 

217 When the two halves are put together, the four sectors corresponding to the 

218 ±15° region around the beam plane cannot be instrumented - see Fig. [31 This 

219 is not a flaw in the design as most of the machine background is located in this 

220 area: should they have been included, these most-lateral sectors would have 

221 been copiously hit by spurious particles. Having sectors 30°-wide in azimuth 

222 is a trade-off between the expected flexibility of the device (in particular one 

223 wants to be able to exclude temporarily a well-defined part of the acceptance 

224 in case of shielding modifications or changes in the background conditions) 

225 and the final number of channels to be instrumented. 

226 Another major constraint of the calorimeter design is the presence around the 

227 IP of the SIDDHARTA detector and of its shielding. Lead bricks are installed 

228 to stop as much radiation as possible before the particles enter the active part 

229 of the experiment. During the initial optimization phase, (February-August 

230 2008) the shielding has been modified to improve the Kaon to background 

231 ratio. The vertical position of this detector has also been tuned (from 5 to 

232 15 cm of the IP) in the meantime. The final shape has been installed in 

233 September 2008. 

234 Each calorimeter sector is a sandwich of twelve trapezoidal 1 cm thick scin- 

235 tillator tiles, wrapped with Tyvek™ [18] paper, alternated with eleven lead 

236 plates of variable thickness. Starting from the face close to the IP, the eight 
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Fig. 3. Drawing representing one arm of the Bhabha calorimeter. One can clearly 
see its structure and its granularity; among the twelve azimuthal sectors, ten are 
instrumented whereas the two around the beam plane are not. 

237 first plates are 0.5 cm thick while the last three are 1 cm thick. Hence the total 

238 thickness of the sector is 19 cm. This design is a compromise between the need 

239 of having a good longitudinal containment of the 510 MeV charged particle 

240 showers (the total depth corresponds to about 12.5 X ) and the requirement 

241 of having a detector not exceeding the permanent quadrupole length - as some 

242 mechanical structures are located right behind these magnets. Moreover, the 

243 further the calorimeter from the IP, the lower the polar angle up to which it is 

244 sensitive. Finally, the free space between the calorimeter and SIDDHARTA is 

245 needed to accommodate the GEM tracker which is described in the following. 

246 The two-hundred and forty scintillator tiles have been produced with injection- 

247 molded technique in IHEP, Protvino [19]. Each tile has three, 2 mm deep, 

248 radial grooves on one face, (one in the middle and two 1 cm away from the 

249 edge of the tile) inside which wavelength shifting (WLS) fibers [20] of 1 mm 

250 diameter are placed - see Fig. H] for details. The thirty-six WLS fibers of a 
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Fig. 4. Picture of calorimeter tiles with the WLS fibers already glued on their radial 
grooves. 

251 given sector are plugged to an optical adapter to fit the photo-cathode of a 

252 Photonis-Philips XP 2262B photomultiplier (PMT) [21] readout by a KLOE- 

253 based data acquisition system described in Section HI 

254 The design sector energy resolution of 35% at 510 MeV is accurate enough 

255 to rely on a threshold on the reconstructed energy to select candidate Bhabha 

256 events. After assembly, the four modules have been tested and calibrated at 

257 the DA$NE Beam Test Facility [22]. Several position and energy (linearity 

258 and resolution) scans have been performed with beams of energy 141, 189, 

259 283, 377 and 471 MeV respectively. Additional pictures of the calorimeter 

260 building process can be seen in Fig. 



26i 3.2 GEM trackers 



262 Rings of GEM detectors [23] aiming at identifying charged particles from 

263 Bhabha events have been designed to be installed in front of each calorimeter, 

264 at a distance of 18.5 cm from the IP. Both trackers are divided into two verti- 
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Fig. 5. Some pictures of a calorimeter module (five active sectors) assembly which 
took place in December 2001 and January 2008 at Frascati. The left picture shows 
the tiles wrapped in Tyvek™ sheets to improve the light collection efficiency; the WLS 
fibers, already bundled together and waiting for being connected to the PMT, are also 
apparent. In the right picture, the module has been put in its holding structure and 
is in the final building steps before being moved to the DA&NE IP. The side view 
allows one to see the twelve scintillator planes such as the eleven lead layers. 




Fig. 6. Assembling a GEM tracker. 

265 cal half moon-shaped units which surround the beam pipe. The top (bottom) 

266 half covers azimuthal angles between 14° and 166° (194° and 346°) respec- 

267 tively. Each of the four GEM trackers is segmented into thirty-two pads: four 
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268 6.5 mm- wide rings contain eight cells covering 19° in azimuth each; the ring 

269 boundaries are located at 68, 74.5, 81, 87.5 and 94 mm respectively from the 

270 beam axis. The trackers are made of three GEMs, thin (50 jum) kapton foils 

271 sandwiched between two copper clads and perforated by a dense set of holes 

272 (70 /im diameter, 140 fim pitch). They are glued on a G10 frame and assem- 

273 bled between the cathode and the anode forming four gaps: from the IP to the 

274 calorimeter, the drift gap (3 mm), two transfer gaps (1 and 2 mm) and the 

275 induction gap of 1 mm. All the volume between the anode and the cathode is 

276 filled with a gas mixture of Ar - C0 2 - CF 4 (45%-15%-40%) [23! . 

277 When a charged particle crosses the drift gap it generates 10-15 primary elec- 

278 trons which drift thanks to the 3 kV/cm electric field before being multiplied 

279 by the GEM foils. As a high potential difference (about 400 kV) is applied 

280 between the copper sides, the holes act as amplifiers and the gain of each 

281 layer is about 20 - hence a factor 8 000 in total. The electron cluster coming 

282 out from the last GEM foil induces a signal on the anode which is amplified 

283 and discriminated by a Carioca Chip [25] plugged on the back side of the pad 

284 PCB. The two front-end electronics cards designed and developed in Frascati 

285 are supplied through a mother board that houses two output connectors with 

286 1 6 LVDS channels each and one input connector with low voltage supplies and 

287 two thresholds channels. High voltage is supplied to the four GEM trackers 

288 with a new system made of four active dividers [26] . 

289 As described above, Fig. [6] shows the different layers of a GEM tracker (left 

290 picture) such as different steps of the assembly process (the four pictures on 

291 the right). Fig. [7] presents one drawing of a calorimeter side with two GEM 

292 trackers in front of it (left plot) and an actual picture (right plot) of the 

293 IP taken at the beginning of the 2008 DA$NE run, when the trackers were 
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294 installed in their nominal location. 

Electron Positron 
monitor Monitor 




Fig. 7. Left: this drawing shows one side of the calorimeter with the GEM trackers 
visible - the rings around the beam pipe in front of the calorimeter structure. Right: 
top view of the DA&NE IP with the SIDDHARTA detector (the three oval-shaped 
structures in the center) and the GEM trackers (the circular shapes in front of the 
calorimeter visible on either side) installed in front of the Bhabha calorimeter. The 
numbering convention for the GEM trackers is described on this picture: trackers 1 
and 2 monitor the electron beam while trackers 3 and ^ are sensitive to positrons; 
odd (even) modules are above (below) the beam line. 



295 3.3 Gamma monitors 



296 Two gamma monitors shown in Fig. [8] are located 170 cm away from the IP 

297 on either side and collect photons radiated by the electron or positron beams. 

298 These detectors are made of four PbWCU crystals which have a rectangular 

299 prism shape: a square section of 30 x 30 mm 2 and a height of 110 mm. The 

300 crystals are maintained together along their longer sides in order to form a 

301 structure 30 mm-wide in front of the photon beam and 120 mm-thick (about 

302 13 Xq). Each crystal is readout by a Hamamatsu R7600 compact photomul- 

303 tiplier [27] which is in direct contact with the top face of the crystal. This 
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304 design helps minimizing the background from particles hitting directly the 

305 PMT. Because of the boost introduced by the IP crossing angle, the photon 

306 trajectories are shifted toward the inner side (along the — x direction) of the 

307 machine. To account for this deviation, the gamma monitors are located at 

308 x = —5 cm from the beam pipe and rotated by 4° in the horizontal plane with 

309 respect to the beam axis. 

310 Thanks to their high logging rates, these detectors provide quick estimates of 
3n the luminosity and background variations. These quantities are then used for 

312 the machine real-time optimization. They cannot easily provide an absolute 

313 measurement of the luminosity as their counts would have to be corrected by 

314 a factor which is continuously changing as the beam currents and the machine 

315 conditions evolve. However, on a short time scale and as relative luminosity 

316 monitors, those counters have been found extremely useful. 

317 3.4 IP shielding 

318 During the course of the 2008-2009 SIDDHARTA runs, ad-hoc lead shieldings 

319 have been constructed around the IP to protect this sensitive detector against 

320 low energy background photons. Three different types of shielding have been 

321 used, called "Soyuz", "Sputnik" and "Mir" respectively. The Soyuz defines 

322 the low angle edge of the calorimeter acceptance by discarding tracks which 

323 are kicked back in the fiducial volume after showering in the QD0, a process 

324 potentially hard to simulate correctly. The Sputnik is meant to absorb the 

325 background generated by the beam bending inside the QD0 magnets; it sur- 

326 rounds the Soyuz. Finally, the Mir shield completes the shielding setup by 

327 closely surrounding the calorimeter. Fig. [9] shows these different components 
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Fig. 8. The gamma monitors are installed 170 cm away from the IP (top picture, 
taken in the DA&NE hall). The window opened through the lead shielding is aligned 
with the trajectory of the photons produced in radiative Bhabha events. Each detector 
is made of four PbWO^ crystals, wrapped in Tyvek™ foils (bottom left picture) . They 
are put together into a PVC box with the four PMTs whose entry windows are in 
direct contact with the top face of the crystals (bottom right picture). 

328 on a side view of the simulated IP - see Section [6] for details. 

329 The shape and the position of these elements are dictated by the existing 

330 machine optics and have been slightly changing over time. All in all, they 

331 result in a partial decrease of the luminometer acceptance. More importantly, 
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Fig. 9. GEANT3 picture showing a side view of the DA&NE IP shielding surrounded 
by the Bhabha calorimeter. The small green components closer to the beam pipe 
form the Soyuz; the red triangles are the Sputnik while the black rectangular volume 
corresponds to the Mir shielding. 

332 they have a significant impact on the GEM detector location and thus on its 

333 performances, as described in Section 15721 

334 4 Data acquisition and trigger 

335 Most of the front-end (FE) electronics and data acquisition (DAQ) compo- 

336 nents, as well as the PMT HV supplies, are taken from the KLOE experi- 

337 ment [28]. The signals coming from the calorimeter are amplified, inverted and 

338 delivered to a 3-stage splitter (SDS board). The first stage is used for timing 

339 measurements and consists of a constant fraction discriminator which deliv- 
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340 ers a current signal to the FE KLOE TDCs [29] whose resolution is 1.04 ns. 



341 Then, the energy measurement stage includes a 3-pole Bessel filter followed 

342 by a KLOE charge ADC (0.25 pC resolution); the filter is needed to match 

343 the PMT signal bandwidth with the ADC input stage. Finally, the third stage 

344 sums all signals coming from the same calorimeter module and this informa- 

345 tion is used as input to the trigger logic. A candidate Bhabha event is defined 

346 as two high-enough energy deposits in back-to-back modules which are com- 

347 patible in time - see Section [5] for details. 



348 The KLOE FE modules are build around a custom bus, the AUX-bus |30j . 



349 which allows a fast readout of the input data as it only uses the VME standard 

350 for initialization. The AUX interface builds sub-events crate by crate ensuring 

351 that the data are properly synchronized. A trigger-driven readout controller 

352 delivers these frames to the control manager board via another custom bus (the 

353 C-bus) which interconnects different crates. The event data are then stored 

354 into a FIFO memory where they are read using VME 32 bits block transfer 

355 mode [31] by a MVME6100 processor running KLOE protocols. This CPU 

356 sends them to a PC computer via a gigabit Ethernet connection where they 

357 are finally written to disk. A simple graphical user interface (GUI) developed 

358 in JAVA performs the run control functions. The entire acquisition chain is 

359 installed in the collider hall, close to the DA$NE IP. Fig. [10] summarizes the 

360 calorimeter data acquisition and trigger systems. 

361 The gamma monitor signals are split after the PMTs: one half is sent to the 

362 charge ADC of the KLOE data acquisition system (and not used further down 

363 for the measurements), while the other goes to an analog mixer. The analog 

364 sum of the four crystals is then discriminated and the counts are read by the 

365 DA$NE control system via a VME scaler readout module. Finally, the GEM 
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Fig. 10. Overview of the Bhabha calorimeter data acquisition and trigger systems. 
This detector is made of four modules (labeled Ml to M4 ) which are split into five 
azimuthal sectors each. A candidate Bhabha event is a time- compatible coincidence 
between two high-enough energy deposits in back-to-back modules: "(Ml AND M4) 
OR (M2 AND MS)". 

366 signal is triggered by the calorimeter; the data coming from the GEM detectors 

367 are discriminated, delivered to FE TDCs and readout via an AUX-bus. This 

368 scheme guarantees synchronization between calorimeter and GEMs. 

369 5 Performances 

370 5.1 The Bhabha calorimeter 

371 Potential Bhabha elastic scattering events e + e~ — > e + e~ are identified at the 

372 trigger level by requiring two coincident energy deposits above a given thresh- 

373 old in back-to-back modules. The locations of the two sectors with the highest 

374 energy deposits are indeed correlated, as shown on Fig. [TTJ Taking as example 
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Fig. 11. Given the numbering convention of the calorimeter sectors (see left plot), 
the position of the two calorimeter sectors with the highest energy deposits can be 
compared event by event. The result (2D-histogram on the right showing the sector 
number with the second highest energy deposit vs. the sector number with the highest 
deposit) shows that these two sectors are most of the time back-to-back, as expected 
from Bhabha events - note that the z-scale is logarithmic. 

375 the opposite sectors 19 and 4, Fig. [12] shows the ADC distributions of each 

376 sector independently, such as the correlation between these two quantities. 

377 Energetic events corresponding to actual Bhabha events are common to both 

378 sectors. 



379 Summing all the sector contributions from a given calorimeter module, Fig. 

380 shows the distribution of the ADC counts for triggered events in the four mod- 

381 ules. Fitting these distributions by a Gaussian, one gets their energy resolu- 

382 tions. 



a(E) 



E 



with E = 510 MeV 



(3) 



384 The results are 22%, 19%, 20% and 20% for the modules Ml to M4 respec- 

385 tively: the design performances are fulfilled in the whole calorimeter. 
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Sector 4 ADC charge 



Fig. 12. Top left and bottom right plots: ADC distributions of the back-to-back 
calorimeter sectors 19 and 4 respectively. The bumps around 1000 counts correspond 
to the true Bhabha events while events for which these sectors did not trigger are 
associated with smaller energy deposits. Top right plot: a 2D histogram showing the 
correlation between the two sector ADC charges; high energy deposits in one sector 
are associated with deposits of similar magnitude in the other sector, the expected 
signature of real Bhabha events. 
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Fig. 13. Distribution of the ADC counts in the four calorimeter modules from events 
in which they are triggering. Fitting these curves provides the corresponding energy 
resolutions which are around 20% at the DA&NE energy (510 MeV); the calorimeter 
design performance is experimentally achieved in this representative high- statistics 
run. 

386 Accidental coincidences involving background particles hitting part of the 
3B7 calorimeter do contribute to the trigger. In order to convert the calorime- 
3B8 ter rate into an absolute luminosity measurement, one needs to identify and 
3B9 remove this noise contamination, whose level is highly dependent on the actual 
390 machine running conditions and cannot be determined a priori. Therefore, a 
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391 filtering procedure has been implemented (first offline and then directly in the 

392 DAQ farm) to improve the accuracy of the online luminosity provided to the 

393 DA$NE control room. The raw trigger rate, R raw , is corrected by a factor 

394 < CF < 1 to give the measured Bhabha rate R me as which is then scaled by 

395 the Monte-Carlo prediction to compute the luminosity - see Eq.([T)). 

396 Rmeas = Rraw X CF (4) 

397 The value of the correction factor CF is adjusted every 3000 events (at most 

398 every few seconds depending on the beam currents) by looking at the distri- 

399 bution of the time difference At between the triggering modules. As shown 

400 in Fig. [14] obtained from a representative high-statistics run, this distribution 

401 has two main components: 

402 • a narrow Gaussian peak centered at At = and corresponding to genuine 

403 Bhabha events (as demonstrated below), 

404 • sitting on top of a flat background coming from random coincidences be- 

405 tween modules and whose width is determined by the duration of the digital 

406 signals building the coincidence (~ 25 ns). 

407 The pattern is similar for both pairs of modules: M1-M4 on the left, M2-M3 

408 on the right. 

409 Two observations allow us to justify the fact that the peak is produced by 

410 real Bhabha events. First, it disappears when the beams are longitudinally 

411 separated - see Fig. [151 Then, the energy deposited in the triggering modules 

412 is significantly larger for these events than for those which originate from 

413 the background. Using the same high-statistics run, Fig. [16] shows the ADC 

414 distributions in the four modules for all events. Three components are visible. 



28 





I I : n i l i o ■ 1 I :~ ? 



-100 -75 -50 -25 



50 75 100 




I Entries 'i4^57 



50 75 100 

l(l)-l{4) [TDC uoiuils] 




-100 -75 -50 -25 25 50 75 100 

1(2)-1{3) [TDC counts) 




1(2)1(3) [TDC t-ounlsl 



Fig. 14. Distributions of the time difference At between back-to-back modules trig- 
gering in coincidence (note the logarithmic y -scale on all histograms). Left column: 
modules Ml and M4; right column: modules M2 and M3. These histograms clearly 
show two components: a flat background from random triggers plus a Gaussian dis- 
tribution centered at At = which corresponds to genuine Bhabha events - see text 
for details. The width of this peak is a few ns (1 TDC count corresponds to 1.04 ns )- 
The background level is about 10% for the M1-M4 trigger and 5% for the M2-M3 
trigger in this example. 

• The white component peaking at corresponds to events for which the other 
pair of modules (M1-M4 for the M2 or M3 histograms; M2-M3 for Ml and 
M4) triggered. As expected, there is none or little energy in the modules 
which did not 'see' a Bhabha decay. 

• The dark-hatched area contains events for which the module triggered while 
the trigger times were not coincident. These events correspond to the flat 
background shown in Fig. [TH 

• Finally, the light-hatched area shows the triggered events which are coinci- 
dent in time. Still comparing this plot with Fig. [Mj these are the genuine 
Bhabha events plus a small background component (note the logarithmic 
y-axis on both figures) one has to subtract to estimate correctly the lumi- 
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Fig. 15. Upper plot: histogram of the trigger time difference between the back- 
to-back modules M2 and M3 accumulated with beams in collisions; lower plot: same 
histogram, but this time for beams out of collisions (longitudinally separated). The 
in-time peak in the timing difference distribution does correspond to luminosity-in- 
duced Bhabha events. 

426 nosity. 

427 Similar conclusions can be drawn from Fig. [17] in which the charges measured 

428 by two opposite modules are summed. Triggering in-time events are mainly 

429 due to genuine Bhabha whereas the two other categories (triggering out-of- 

430 time or not triggering) are background-related. 



To estimate the amount of background under the Bhabha peak, a sideband 
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Fig. 16. ADC distributions in the four calorimeter modules for the high-statis- 
tics run used to produce most of the plots in this section: all histograms show three 
components with similar patterns. Events for which that particular module did not 
trigger belong to the white area: little energy is deposited in these sectors, as expected 
from background fluctuations. Hence, these events are not real Bhabha decays for 
which one of the two particles is not detected. The hatched regions contain the other 
events for which the module triggered. Events for which the energy deposit is not 
coincident in time with the one measured in the opposite module (pure background 
case) are in the darker area. The deposited energy is lower than for the in-coinci- 
dence events (high purity Bhabha sample) which are in the light-hatched area. These 
plots and those in Fig. \15\ justify the timing criterion used to separate signal from 
background. 
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Fig. 17. Distributions of the sum of the energies deposited in opposite modules 
(M1-M4 and M2-M3) in the high-statistics run taken as example in this section. The 
three categories of events - not triggering (white area), triggering but coincidence 
out of time (dark hashed area) and triggering in-time (light hashed area) - are clearly 
visible and confirm that the Bhabha events fall in the third category. 



432 region containing only random coincidences is defined - see the hatched area 

433 in the bottom plots of Fig. [TH Counting the number Ni n of events in the At 

434 range [-6;6] ADC count (corresponding to about ±3<r of the Gaussian Bhabha 

435 peak, 1 count is 1.04 ns) and the number N s b of events in the sideband region 

436 (of equal width: 12 ADC counts), one has: 



CF = 1.0 - ^ (5) 

Na 
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438 This simple procedure provides an accurate real-time measurement of the 

439 absolute luminosity which is among the main feedback signals used by the 

440 DA$NE operators. Fig. [TS] displays about 35 minutes of data taken on May 

441 24 th 2008 to show the improvement brought by this background-subtraction 

442 algorithm. The top plot displays the beam currents, the middle one the uncor- 

443 rected L uncorr and corrected L corr (background-subtracted) luminosities and 

444 the bottom one the ratio L corr / L uncorr . Apart for a few marginal points, the 

445 uncorrected luminosity is always greater than the corrected one, by 10-20% in 

446 this particular example. The difference is even more striking between minutes 

447 26 and 29, a time period during which the beams are longitudinally sepa- 

448 rated. If L corr goes immediately to as expected, L uncorr remains sensitive to 

449 a background component whose decay time is on the timescale of a minute. 

450 5.2 The GEM trackers 

451 Unfortunately the GEM trackers could only be kept in their nominal position 

452 (18.5 cm from the IP) until April 2008 when the SIDDHARTA lead shielding 

453 had to be extended to protect efficiently this detector from machine back- 

454 ground. In that initial period, the GEM trackers were left out of the DAQ 

455 due to a conversion problem between the 128 LVDS Carioca channels and 

456 the ECL inputs of the KLOE TDCs. However, the four trackers were used as 

457 background monitors sending data every second to the DA$NE control room. 

458 This information was available all the time, in particular during the injection 

459 phases. Fig. [T9l shows a snapshot of the DA$NE online stripcharts monitoring 

460 the background levels measured in the GEM trackers. Pad-to-pad coincidences 

461 (64 opposite couples in OR) were also studied to get another estimate of the 
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Fig. 18. Display of about 35 minutes of data taken on May 24 2008. From top to 
bottom: beam currents, uncorrected L uncorr and corrected L corr luminosities, ratio 
L corr /L uncorr . For this particular period of running, the background level is between 
10-20%. Moreover, when the beams are longitudinally separated (between minutes 
26 and 29) L 

corr goss immediately to (as expected) whereas L uncorr does not. 

462 Bhabha rate. Yet, accidental coincidences were found to exceed significantly 

463 the geniune Bhabha rate and no significant variation was observed when the 

464 beams were put out of collision. 

465 In July 2008 the GEM trackers were installed again for a few days, but only 

466 1 cm away from the IP in front of the lead shielding. As they were finally 

467 included in the DAQ, data were acquired when the calorimeter was triggering. 

468 Fig. [20] shows the efficiency of two of the four GEM modules as a function of 

469 the offline calorimeter ADC threshold; the plateau at 40% is due to the shift 

470 of the tracker ^-positions: as the modules are closer to the IP, a large fraction 
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Fig. 19. These plots show how the GEM trackers were used as background monitors 
at the beginning of the upgraded DA&NE commissioning - the rates presented in this 
example are typical of stable beam conditions. Each line of two plots corresponds to 
one of the four GEM trackers (recall that modules 1 and 2 are sensitive to the 
electron beam, 3 and 4 to positrons). On the left column, rates (in Hz) measured in 
each of the 32 tracker cells (a GEM module has an 8-fold segmentation in azimuth 
and a J^-fold segmentation in radius, see Section \3.2\) are shown: the lighter the 
color, the higher the rate. The background is typically around 1 kHz per pad. On 
the right column, the time evolutions of the tracker rates are displayed. Updated at 
1 Hz, these plots (whose axis boundaries can be setup independently, which explains 
why the settings are slightly different on the four snapshots) have been used in real 
time by the DAQNE operators to monitor the machine conditions. At the time this 
picture was taken, the electron beam background in the GEM was about 28 kHz while 
it was around 48 kHz for the positrons. More background was seen above (below) the 
beam line on the electron (positron) side. 
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472 
473 



of Bhabha tracks hitting the calorimeter pass through the GEM central hole 
and are thus not detected. The average hit cluster size was 1.2 pad; Fig. [21] 
displays the 9 and correlations of the hits in the four GEM detectors. 



474 When the nominal SIDDHARTA data taking run started, the GEM trackers 

475 were definitely removed from the IR. 



200 400 



Fig. 20. Efficiencies (in %) of the two GEM trackers located above the beam line 
(#1 is sensitive to electrons, #3 to positrons) versus the threshold applied offline on 
the calorimeter ADC output: at most 40-45% of the calorimeter tracks are seen in 
the GEMs. This is explained by the fact that the trackers were closer to the IP than 
foreseen: a significant fraction of the tracks hitting the calorimeter passed through 
the GEM center hole and were not seen. The effect is enhanced by the strong polar 
angle dependance of the Bhabha cross-section. 



476 5.3 The gamma monitors 



477 The gamma monitors are mainly used to optimize the luminosity in real time 

478 as they are very sensitive to sudden changes in machine conditions. Thanks 

479 to the high radiative Bhabha cross section at low angle, the photon rate typ- 

480 ically reaches several tens of kilohertz for a luminosity around 10 32 cm~ 2 s . 



36 



bhobho dtetrlbullon 



I ENTRIES an" 



10 12 14 16 




Fig. 21. Example of data acquired by the GEM trackers to check the collinearity 
of Bhabha events triggered by the calorimeter. The top plot shows the azimuthal 
correlation of the charged tracks - each pair of GEM trackers contains 16 (^-sectors, 
8 above the beamline and 8 below, whose numbering is chosen so that opposite sectors 
have the same label. The bottom plot shows the anticorrelation between the track 
polar angles. Combining both information shows that the Bhabha events are back to 
back, as expected. 

48i Therefore, these detectors can provide measurements with negligible statistical 
4B2 fluctuations every second. This is illustrated by Fig. [22] which shows the time 
4B3 evolution of the gamma monitor rates (top plot: e + side, red histogram; mid- 
434 die plot: e~ side, blue histogram) and of the background-subtracted Bhabha 
4B5 luminosity (lower histogram; in green) during a machine test in April 2008. At 
4B6 a certain time clearly visible on the charts, the beams are put out of collision 
4B7 by a 180° RF phase shift; this new condition is immediately detected by the 
4B8 gamma monitors whereas there is some latency in the Bhabha rate. Sharp 

489 variations are also visible when the beams go back in collision. The smooth 

490 decrease of rates over the whole time range, visible in all three plots, corre- 
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Fig. 22. Time evolution of the gamma monitor (top and middle charts, e + side in 
red and e~ side in blue) and Bhabha background- subtracted (bottom chart in green) 
rates recorded during a machine test performed on April 22, 2008. All rates drop 
drastically when the two beams are separated by a 180° shift of the RF phase although 
the gamma monitors react quicker than the calorimeter. Similar conclusions can be 
drawn when the two beams collide again. 

491 sponds to a period of coasting for both DA$NE beams. One can also note 

492 that the residual background levels are very low in the absence of collision. 

493 The quick response of the gamma monitors and their low background con- 

494 tamination levels allow them to be used for precision measurements of the 

495 transverse beam size. The scan is carried out by shifting horizontally and ver- 
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Fig. 23. DA&NE control room monitor showing the result of a beam vertical scan 
(moving the electron beam while keeping the positron one fixed) performed in April 
2008. Fitting the data coming from the gamma monitor allows one to estimate the 
quadratic sum of the two beam vertical sizes. The result, as displayed in the top box 
on the right column, is Y>™ eas = 9.34 fim for this particular scan. 



tically the beam trajectories. The rate variations measured during such test 
give access to the bunch profile. In order to be successful, this procedure must 
be very fast to avoid the result to be biased by the contribution of other ef- 
fects (e.g. loss of luminosity when the beams are coasting, an effect which 
is relevant in the DA$NE ring where typical beam lifetimes are lower than 
1000 s). Fig. [231 shows a typical vertical scan performed in April 2008. In this 
example, the measured quadratic sum of the two beam vertical sizes is found 
to be Y^ ea& = 9.34 fim. To match the upgraded DA$NE running conditions, 
computations [32] have shown that a corrective factor of 0.88 needs to be ap- 
plied to get the real value of S r Assuming that the two beams are Gaussian 
and of equal sizes, one finally gets o y = 5.8 /^m for this particular scan. 
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507 6 Monte-Carlo simulation 

508 A Monte-Carlo simulation of the luminometers is needed to validate the (ra- 

509 diative) Bhabha selection algorithms, to find the correspondence between the 

510 measured event rates and the actual machine luminosity and to study the 

511 impact of the background. For this purpose, a package has been developed to 

512 simulate the detector response in the GEANT3 [33] framework. Input events 

513 are either Bhabha generated with the BHWIDE [34J Monte-Carlo generator 

514 or simulated particles leaving the machine nominal orbit close to the IP |14j . 

515 After a detailed description of the simulation code (from the generation to 

516 the reconstruction of the events in the virtual detector), the main simulation 

517 results are presented. 

sis 6.1 Event Generation 

519 The simulation is a two step-process. Events are first created at the genera- 

520 tor level using the BHWIDE code: both non-radiative and radiative Bhabha 

521 events (without limitation on the number of radiated photons) can be simu- 

522 lated. The BABAYAGA package |35j has provided cross-checks of the BH- 

523 WIDE performances as this software had been used in the past by the 

524 KLOE collaboration for a precise determination of the luminosity delivered 

525 by DA$NE between 2000 and 2007 [36]. Both generators have been found 

526 in agreement. The four- momenta of the generated particles (assumed to be 

527 produced at the IP), as well as those of the potentially initial-state-radiated 

528 photons are stored in a PAW [37] ntuple which is written to disk in order to 

529 be used as input to the GEANT simulation. 
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530 Like for the real machine (see Fig. [2]), the generated positrons are pointing 

531 toward the direction of the +z axis while the electrons fly in the opposite 

532 direction. Convenient software switches allow one to select particular polar 

533 angle 9 ranges for the charged tracks (whose distribution scales like 1/8 3 in 

534 the CM frame). The horizontal boost due to the 25 mrad crossing angle be- 

535 tween the two incoming beams and pointing toward the — x direction is also 

536 simulated. As shown on Fig. [2H its effect is quite significant: the polar and 

537 azimuthal correlations between the Bhabha positrons and electrons are mod- 

538 ified by the non-zero crossing angle. Although it may look small, it impacts 

539 significantly the correlation between the particle trajectories which makes the 

540 understanding of its consequences mandatory to ensure accurate luminosity 

541 measurements. 

542 Particles leaving the nominal orbit close to the IP are expected to be the main 

543 source of background for the luminosity measurements. They are simulated 

544 with a code developed for the DA$NE experiment [H] which parameterizes 

545 the particle coordinates as a function of the path length (s) along the design 

546 trajectory. Touschek events are generated continuously all over the ring, av- 

547 eraging the probability density function every three machine elements. Then, 

548 particles are tracked according to the optics over many turns or until they 

549 are lost. The actual beam pipe aperture is taken into account; nonlinear ele- 

550 ments, sextupoles and octupoles, are also included in the tracking, allowing an 

551 intrinsic dynamical aperture calculation. As this simulation does not include 

552 the experimental setup, a dedicated procedure - see next section for details - 

553 needs to be applied prior to the GEANT simulation, in order to go from the 

554 beam path frame to the DA$NE coordinate system. 
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Fig. 24. Simulated Bhabha events without (top row) and with (bottom row) the 
DA&NE horizontal crossing angle. CM generation is done with BHWIDE before 
the track parameters are converted to the laboratory frame described previously. 
The left column compares histograms of the sum of the two track polar angles: 
£0 = # P os + #eW the right column shows the corresponding histograms for the vari- 
able A<£ = 1 P os — 0cic I describing the azimuthal correlation between the electron and 
positron tracks. The S© and A<£ distributions vary in shape with the boost and the 
mean value of the latter is also modified - beam trajectories are pointing toward the 
—x direction. The non-zero widths of the top (no boost) distributions are due to the 
final state radiated photons. 

555 6.2 Detector Simulation 

556 The ntuple produced by BHWIDE is used as input of a FORTRAN-based 

557 GEANT3 code which simulates the region around the IP relevant for the 

558 luminosity measurements: from the e + gamma monitor to the e~ one (z = 
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559 ±170 cm). The generated particles are positioned at the IP from where 

560 GEANT propagates them until they leave the simulated mother volume, decay, 

561 or get absorbed. Along the way, GEANT computes the effect of their interac- 

562 tions in the various materials they cross. Information relevant for the study of 

563 reconstruction algorithms or for the detector design optimization (such as the 

564 true impacts in the GEMs and colorimeter modules, or the number of photons 

565 detected in each calorimeter sector or in the gamma monitors) is stored in an 
see output PAW ntuple for further analysis. 

567 Fig. [25] gives an overview of the experimental setup as implemented in the 

568 GEANT simulation. The left part shows the setup from above the beam 

569 plane - the SIDDHARTA volume has been removed to make the drawing 

570 clearer. The main magnified elements are the Bhabha calorimeter modules, the 

571 gamma monitors, the QDO permanent quadrupole magnets and the vacuum 

572 chambers whose geometry is particularly complicated in the region where the 

573 beam pipe, common to both beams at the IP, gets split (Y-tubes) to allow par- 

574 tides of opposite charges to be separated in the rest of the DA$NE ring. The 

575 top right drawing shows details of the Soyuz, Sputnik and Mir lead shielding 

576 elements (see Section |3~I4"|) ; for a better readability, the Bhabha calorimeter has 

577 been removed from this picture. The big shieldings protecting the calorimeters 

578 from backward-orientated particle showers are also visible. Finally, the bot- 

579 torn right picture is a cut in the y-z plane which explains how crucial it is to 

580 simulate these shields accurately. The red tracks are Bhabha generated with 

581 the BHWIDE package whose density reflects the l/# 3 -dependence of the cross- 

582 section. Dramatic rate variations are visible around the shield edges, which 

583 means that a small error in their position would trigger large acceptance effects 

584 as the fiducial Bhabha volume would be incorrect. 
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Fig. 25. Overview of the experimental setup as implemented in the GEANT simu- 
lation; see text for details. 
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585 The various luminometer elements, the machine components and shielding 

5B6 relevant for the luminosity measurements are simulated, with the proper di- 

5B7 mensions and materials - see Fig. 

5B8 • The scintillator tiles, wrapped in their Tyvek™ skin, and the lead absorber 

5B9 layers of the calorimeter. Losses of light in the tiles and in the optical fibers 

590 connected to the PMTs are also taken into account, such as the PMT quan- 

591 turn efficiencies. 

592 • The GEMs, with both their azimuthal and radial segmentations; the three 

593 kapton foils and their copper coating, the four longitudinal layers of gas 

594 (one drift and three transfer volumes) as well as the copper box containing 

595 the whole tracker. 

596 • The gamma monitor, made of four PbW04 crystals wrapped in a Tyvek™ 

597 skin and put together in a PVC box. 

598 • The various sections of the aluminum beam-pipe, including the splitters 

599 (the green volume downstream of the calorimeter on the top-right plot) and 
eoo the three cylindrical sections. 

601 • the QF1 (the cylinders upstream of the splitters) and QDO quadrupoles 

602 (invisible as surrounded by the tile calorimeters), whose magnetic fields are 

603 also simulated. 

604 • The Soyuz, Sputnik and Mir shields. 

605 Like for the Bhabha event generation, the input parameters are readout from 

606 a datacard and several software switches have been implemented to allow the 

607 user to customize the simulation: value of the horizontal beam crossing angle 
60s (defining the size of the boost), choice of the hardware pieces to be included in 
609 the IR and of the particles (electrons, positrons or photons) to be propagated 
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Fig. 26. Distribution of particles lost in the IR versus path length s (top plot) and 
the corresponding trajectories (bottom plot). The IP is at s = and the collimators 
are inserted; the blue and red rectangles show the position of the QDO and QF1 
quadrupoles respectively. 

6io in GEANT, etc. 



en As mentioned in the previous section, the simulation of the particles lost from 

612 the beams - see Fig. [26]- does not include the experimental setup around the 

en IP. Hence, one needs to process the output of this code (positions and momenta 

614 of the particles when they leave the nominal orbit) before using it as input to 

615 the GEANT simulation. The first step of this procedure consists in converting 

616 the particle parameters from a frame linked to the beam design trajectory to 

617 the IP frame defined in section El Knowing the path length s, the particle 
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618 transverse coordinates (x s ,y s ), their first derivatives with respect to s (x' s ,y' s ) 

aw and the energy deviation with respect to the nominal value of 510 MeV, one 

620 can compute the particle position and momentum in the (x, y, z) system of 

621 coordinates. 

622 Then, these particles are propagated backward in time until they are inside 

623 the beam pipe. For this purpose, one flips the true particle momentum and 

624 uses an 'inverted' GEANT simulation of the IR where the vacuum beam pipe 

625 is now lead and all other materials vacuum. With this setup, the first GEANT 

626 hit associated with a non-zero energy loss can only occur inside the beam pipe 

627 which provides the true position of the particle before it would leave it. Its true 

628 momentum is then defined as the opposite of the value computed by GEANT 

629 at the last hit with null energy loss. These parameters are finally stored on disk 

630 in another ntuple which can then be used as input of the 'regular' GEANT 

631 simulation which is finally used to simulate the Touschek background rates 

632 in the gamma monitors. As can be seen on Figs. [8] and [251 these detectors 

633 are surrounded by lead bricks. This shield allows them to fulfill their primary 

634 task: a clean measurement of the photon emitted by radiative Bhabha events, 

635 in spite of the large beam backgrounds. Consequently, the acceptance of the 

636 gamma monitors is limited to a circular window, dug in the lead and aligned 

637 with the typical line of flight of a radiative Bhabha photon, whose diameter 

638 does not exceed 1 cm. Because of this reduced acceptance, Touschek particles 

639 are not expected to directly hit the gamma monitor. They have to be measured 

640 indirectly, via the secondary particles produced in showers caused by their 

641 passage through the QF1 magnet or the beam wall. A significant fraction 

642 of these secondaries have their energy lower than 10 keV. This is the limit 

643 below which GEANT is no longer reliable in its description of the particle- 
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644 material interactions. This is why the simulation of the Touschek background 

645 measurement is based on the GEANT4 [38] package, which has been improved 

646 w.r.t GEANT3 to better deal with the low energy regime. 

647 6.3 Reconstruction 

648 Let's now focus on the reconstruction of events in the Bhabha calorimeter. 

649 GEANT3 simulates the showers initiated by the charged particles and com- 

650 putes in particular the amount of energy lost at each step of the tracking while 

651 the particles are crossing the scintillating tiles. The number of recorded pho- 

652 tons is deduced from this quantity using a set of realistic parameters which 

653 are provided as input to the simulation: the scintillator photon yield per unit 

654 energy; the attenuation depth in this material and in the WLS optical fibers; 

655 the PMT quantum gain. The top plot in Fig. [27] shows the distribution of 

656 the number of reconstructed photons per event in a given module, summing 

657 up the 'most hit' sector and its two closest neighbors. As the average number 
ess of photons correspond to the nominal beam energy of 510 MeV, the energy 

659 resolution appears to be 18%/y/E. 

660 Fig. [27J also shows the dependence of the number of collected photons on the 
eel true track angular coordinates (8 and <ft) at the calorimeter. The </> depen- 

662 dence (bottom left plot) comes from the calorimeter segmentation and from 

663 the correlation between the azimuthal track position and the distance be- 

664 tween the hits in the tile and the fibers collecting the light in the material. 

665 This correlation also explains the 9 dependence visible on the bottom right 
eee plot. The amplitude of this effect is obviously function of the properties of 
667 various components: scintillator, WLS fiber and PMT. These topics are dis- 
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Fig. 27. The top plot shows the distribution of the total number of photons collected 
by the PMTs in the "most hit" calorimeter sector and its two closest neighbors. The 
histogram has a Gaussian-like shape whose resolution is about 18%, in agreement 
with the calorimeter design and the actual data. The two bottom plots show how this 
number depends on the true azimuthal <f> (left side) and polar 9 (right side) track 
angles. On the left plot, the calorimeter segmentation is clearly visible: the tracks 
which deposit the highest energy in average are those which hit the middle of the 
azimuthal sectors. In addition, the right profile histogram shows that the deposited 
energy decreases with the polar angle as the tracks get less contained. 



ees cussed in details in Section [7] which is dedicated to the validation studies of 
669 the simulation and where a good agreement between data and Monte-Carlo 
e/o distributions is shown. 
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671 6.4 Calorimeter thresholds 



672 Another challenging part of the Monte-Carlo simulation consists in choosing 

673 a proper set of thresholds to control the individual module triggers. As the 

674 reconstructed energy is directly proportional to the number of photoelectrons 

675 at the photomultiplier, the trigger thresholds are applied at this level. This 

676 is not an accurate description of the calorimeter trigger: in the real system 

677 (see Section S]), each module threshold is applied on a pulse coming from the 

678 KLOE SDS board summing up the currents collected by the PMTs in the five 

679 sectors. This electronics chain adds some noise that smears the signal; thus, 
6bo the pulses are not directly proportional to the numbers of photoelectrons. A 
6bi systematic uncertainty is estimated to account for this discrepancy between 
6B2 data and Monte-Carlo. The effect of the smearing is evaluated by varying 
6B3 the thresholds used in the simulations within a range centered on the value 
6B4 determined from the data (see below), and twice as large as the error on this 
6B5 determination. 

686 The Monte-Carlo thresholds are adjusted to those found in data which are 

687 computed using the following procedure. Taking the distribution of the total 
ess energy in a given module (defined as the sum of the pedestal-subtracted ADC 

689 values of the five associated sectors) and a cut value, one can count the number 

690 of events for which that energy exceeds the cut. Then, one can compute the 

691 fraction of those events for which the module actually triggered. The curve 

692 showing the variation of this fraction versus the cut value has a step-like 

693 shape which allows one to extract the threshold. For low (high) cut values, 

694 the fraction is close to (1); the transition between the two states is sharp 

695 and occurs around the threshold. For the reason explained in the previous 
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Fig. 28. Fits of the fractions of triggering events vs. ADC cuts for the four calorime- 
ter modules. The captions show the fit results, as described in the text. 

696 paragraph, the curve is not a pure step function but exhibits some Gaussian 

697 smearing. Therefore, it is fit with a 3-parameter model: 



(c, p,r,a) = p x E 



c — r 



a 



(6) 



699 where 



700 • c is the cut value; 

701 • p is the ratio plateau (in the high cut region) which should be very close to 

702 1; 

703 • r is the fitted threshold; 

704 • a is the smearing parameter, i.e. the width of the Gaussian; 

705 • finally, E is the Gaussian error function. 

706 Fig. [28] show the fit results for the four calorimeter modules, using as input 

707 the data from a representative DA$NE run. The uncertainty in the threshold 
70s determination is mainly related to the background contribution which distorts 

709 the left tail of the fitted function. The systematics effect, estimated by varying 

710 the low boundary of the fit range, is around 5%. 
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711 Finally, the ADC thresholds need to be converted to energy values. As gen- 

712 uine Bhabha events deposit in average 510 MeV in triggering modules, the 

713 peak of the ADC distribution obtained by only keeping events for which the 

714 module triggered provides the required conversion factor, with a systematic 

715 uncertainty of the order of 1% (again background-driven). Table [1] summarizes 

716 the threshold values of the four calorimeter modules, using as ADC-to-energy 

717 conversion factors the fit results shown in Fig. [13j 



Module ADC Threshold value r Energy 

1 301.9 135 MeV 

2 346.7 151 MeV 

3 318.9 146 MeV 

4 311.1 145 MeV 

Table 1 

ADC count and energy thresholds for the four calorimeter modules. 

718 6.5 Expected Bhabha rates 

719 Finally, Table [2] summarizes the predictions of the Monte-Carlo simulation 

720 for the event rate in the Bhabha calorimeter for various shielding configura- 

721 tions. These rates correspond to a benchmark luminosity of 10 32 cm~ 2 s _1 . The 

722 systematics uncertainties which are quoted here are computed in Section [7J 

723 These rates are high enough to allow an online measurement refreshed every 

724 few seconds with a small statistical uncertainty. 
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Soyuz only 


Soyuz + Sputnik 


Final Setup 


Bhabha Calorimeter 


470 ± 65 Hz 


280 ± 40 Hz 


230 ± 30 Hz 



Table 2 

Event rates predicted by the simulation in the Bhabha calorimeter, assuming a 
luminosity of 10 32 cm _2 s _1 . 

725 7 Validation studies 

726 This section deals with the evaluation of the systematic uncertainties for the 

727 online luminosity measurements performed with the Bhabha calorimeter. It 

728 also presents the data-Monte-Carlo comparison studies which have been car- 

729 ried out to assess their validity. 

730 7.1 Systematic uncertainties on the calorimeter-based luminosity measure- 

731 ments 

732 The uncertainty on the Bhabha luminosity measurements stems from a variety 

733 of sources. Most of them are related to inaccuracy or uncertainties in the 

734 Monte-Carlo simulation. 

735 7.1.1 Bhabha generator uncertainties 

736 The systematic uncertainty due to potential inaccuracies in the Monte-Carlo 

737 kinematic variable distributions and differential cross-sections is given by the 

738 variation in the Bhabha rate predicted by the simulation when using the 

739 BABAYAGA event generator instead of BHWIDE. It is found to be small. 
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740 7.1.2 Uncertainty on the interaction region phase space 

741 The distributions of the position of the Bhabha vertices and of their CM 

742 momenta depend on the phase space distribution of the colliding bunches, on 

743 the hourglass effect [39] and on the beam crossing angle. A precise prediction of 

744 the shape of those distributions is not an obvious task, made even more delicate 

745 by the crab waist scheme. Writing a realistic IP phase space simulation would 

746 be very time-consuming. However, this is not necessary as the crab waist 

747 scheme reduces the impact of the hourglass effect by making the length of 

748 the overlap between two colliding bunches of the same order as the vertical 

749 betatron function. Consequently, the hourglass effect is not simulated: the IP 

750 distribution used in the Monte Carlo simulation, is a mere three dimensional 

751 Gaussian, whose standard deviations are taken as the quadratic sum of the 

752 colliding bunch sizes. The longitudinal sizes are reduced to account for the 

753 crossing angle. The uncertainty associated with this approximate description 

754 is estimated in a conservative way by quoting the variation observed in the 

755 Bhabha rate when the length of the longitudinal component of the Gaussian 

756 is set to the value corresponding to a null crossing angle. 

757 7.1.3 Geometrical uncertainties 

758 The location of the various hardware elements (detectors, shields, etc.) defining 

759 or impacting on the Bhabha acceptance is known with a precision of the order 

760 of a pos = 2 mm. We vary these positions by ±cr pos and rerun the Monte-Carlo 

761 simulation using these slightly modified setups. The corresponding variation 

762 in the acceptance is taken as the associated systematic uncertainty. It is the 

763 dominant contribution to this error and its main contribution comes from the 
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764 uncertainties in the SIDDHARTA shield position. 

765 7.1.4 Uncertainty on the simulated energy reconstruction 

766 The procedure followed in the Monte-Carlo simulation to derive the recon- 

767 structed energy from the deposits associated with the GEANT3 raw hits is 

768 described in Section [61 It is the source of several systematic uncertainties 

769 which are reviewed in the following. The photon yield in the scintillator fol- 

770 lows a Poissonian statistics whose parameter depends on the nominal yield 

771 per unit energy which is not perfectly known. Consequently, the resolution on 

772 the reconstructed energy - and thus the simulated acceptance of the trigger - 

773 could be wrong. The energy resolution is also impacted by the noise level in 

774 the PMT amplification. Another contribution to this uncertainty comes from 

775 the absorption of the scintillation photons while they are traveling through the 

776 tiles or along the WLS fibers. For instance, two equivalent energy deposits can 

777 yield different reconstructed energy values if one is closer to a fiber than the 

778 other. This makes the average reconstructed energy dependent on the track 

779 polar and azimuthal angles, as shown in Fig. [27J This dependence broadens 

780 the overall reconstructed energy distribution in a way that has to be correctly 

781 taken into account if one wants to reach an accurate simulation of the energy 

782 resolution. Since the attenuation length of the materials the scintillator tiles 

783 and WLS fibers are made of is not precisely known, it was not the case in the 

784 first place. This problem has been solved thanks to the test beam data: plots 

785 similar to those shown on Fig. [57J have been produced and used as reference 

786 to tune the Monte-Carlo attenuation length. 

787 Finally, in Section 17.41 a quite good agreement between the actual Bhabha 
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7B8 event energy distribution and the corresponding Monte-Carlo prediction is 

7B9 shown. This explains why we eventually quote no systematics associated with 

790 the energy reconstruction. 

791 7.1.5 Uncertainties due to the electronic threshold determination 

792 The treatment of the trigger by the simulation is described in section 16.31 

793 The energy thresholds are varied within ranges determined from the data. 

794 The associated variation of the Bhabha rates is taken as the corresponding 

795 systematic uncertainty. 

796 7.2 Background-related uncertainties 

797 The background treatment is described in Section [5j Two sources of systematic 

798 uncertainties associated with the real-time background subtraction procedure 

799 have been identified. First, the number of events in the sideband of the time 
boo difference distribution between two back-to-back triggering module exhibits 

801 statistical fluctuations which are directly translated into a systematic error on 

802 the corrected event rate. Secondly, the procedure also assumes that the peaking 

803 background is negligible. The validity of this assumption has been verified by 

804 checking that the calorimeter luminosity variations are in agreement with 

805 the expectations when beam currents and sizes vary. Moreover, the genuine 

806 Bhabha timing peak disappears completely when the beams are put out of 

807 collision - see Fig. [151 

80s For completeness, we also run the Monte-Carlo simulation on two-photon 

809 events generated by BHWIDE as their back-to-back topology makes them 
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810 an irreducible background. The low rate predicted by the simulation (about 

en 0.1%) is subtracted from the number of signal events obtained after the online 

812 subtraction procedure. The same number is conservatively quoted as system- 

813 atics. 



814 7. 3 Results 

815 Table [3] summarizes the various contributions to the systematic uncertainties 
aw presented in the previous paragraphs. The total systematic uncertainty on the 
si/ measured luminosity, taken as the quadratic sum of all these contributions, is 
sis of the order of 15% and is dominated by the geometrical uncertainty. 



Source 


Value % 


Bhabha Generation 


2% 


IP Phase Space 


4% 


Geometry and Alignment 


11% 


Threshold Determination 


5% 


Background Treatment 


3% 


Total 


13% 



Table 3 

Systematic uncertainties on the calorimeter luminosity measurements. 
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Fig. 29. Trigger rates per sector for the four calorimeter modules, with the final 
setup. Red points are the Monte-Carlo expectations while black histograms are from 
real data. The rate variations between sectors are due to the boost and to the angular 
coverage of the various shields. 

8i9 7.4 Data - Monte-Carlo comparison 



820 The predictions of the Monte-Carlo simulation for a number of quantities have 

821 been compared with the data to verify that no effect larger than the systematic 

822 errors described in the previous section has been overlooked. 

823 The Bhabha rate per calorimeter sector is shown on Fig. [291 both for data 

824 and Monte-Carlo. This distribution is not flat because of the boost and of the 

825 the various IR shieldings. The Monte-Carlo-data agreement is satisfactory al- 

826 though the various sectors are not equally affected by simulation inaccuracies. 

827 This figure also shows that all the sectors work properly in the real detector 
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Fig. 30. Trigger rates per sector for the four calorimeter modules before the installa- 
tion of the SIDDHARTA final shield. Red points are the Monte- Carlo expectations 
while black histograms real data. The rate variations between sectors are due to the 
boost, and to the angular coverage of the various shields. 

828 (no broken fiber, no module assembly issue, no problems in the PMT connec- 

829 tion quality or isolation from outside light, etc.). It is interesting to notice that 

830 a similar agreement was observed between data and Monte-Carlo during the 

831 period preceding the SIDDHARTA shield installation (see Fig. 130]) . Obtaining 

832 this result for two different setups is an indication that the overall acceptance 

833 is correctly described. 

834 The effects mentioned above, as well as the imperfect knowledge of the attenu- 

835 ation length in the tiles and WLS fibers also affect the energy reconstruction, 

836 which can be tested in detail by studying the resolution on the electron or 

837 positron azimuthal angle, reconstructed via a 'barycenter' method (average of 
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Fig. 31. Resolution on the azimuth angle, reconstructed as the weighted average of 
the <j) positions of the calorimeter sector centers, the weights being the values of 
the energy deposited in each sector - this method is called 'bary center' in the text. 
The resolution is taken as the difference between the electron and positron 4> angles. 
The shape asymmetry seen in both plots (top:Ml-M4 coincidences; bottom: M2-M3) 
is due to the boost caused by the beam crossing angle. The missing Monte- Carlo 
points in the left part of these plots is due to the logarithmic scale and to the limited 
statistic of the Monte- Carlo samples which have been rescaled to the data sample 
statistics to ease the comparison. 

838 the hit azimuth angles using the corresponding energy deposits as weights). 

839 The resolution can be estimated by subtracting the electron reconstructed 
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840 azimuth <p from that of the positron. In the absence of boost, this quantity 

841 only differs from zero because of resolution effects. The distortion due to the 

842 boost makes this variable also useful to check that it is correctly implemented 

843 in the simulation. Fig. [31] compares the azimuthal resolution in Monte-Carlo 

844 with the one extracted from data. The agreement, although not perfect, is 

845 reasonable. Comparing the shape of the Bhabha energy peaks in data and 

846 Monte-Carlo (Fig. [32]) show that the discrepancies observed on Fig. [31] do not 

847 impact enough the energy distribution to bias the trigger acceptance. It also 

848 validates the simplified description of the trigger thresholds in the simulation. 
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Fig. 32. Comparison of the reconstructed energy distributions for data (black his- 
togram) and the Monte- Carlo simulation (red dots) in the four calorimeter modules: 
Ml (top left), M2 (top right), M3 (bottom left) and M4 (bottom right). 
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849 8 Results 



850 The Bhabha calorimeter and the gamma monitors have been designed, built 

851 and operated to monitor the expected improvements of the DA$NE machine 

852 after its 2007 upgrade and to provide a fast and useful feedback to the control 

853 room operators. These goals have been achieved demonstrating that the ac- 

854 celerator operations are successful. The last section of this paper reviews the 

855 main results of the 2008-2009 DA$NE runs which reflect the performances of 
sse both the machine and the associated detectors. 

857 8. 1 Luminosity 

858 Fig. [33] shows a typical two-hours running period at DA$NE. The peak cur- 

859 rents are around 2 Amp for electrons (top plot, blue curve) and 1 Amp for 

860 positrons (top plot, red curve) while the instantaneous corrected luminosity 

861 (middle plot, red dots) peaks between 3.5 and 4.0 10 32 cm" 2 s _1 . As expected, 

862 backgrounds (in particular injection- related) only affect the uncorrected lu- 

863 minosity (middle plot, green curve). Although an high-injection regime is not 

864 compatible with the operations of the background-sensitive SIDDHARTA ex- 

865 periment, this running mode has been tested for few hours; integration in 
see excess of 1 pb _1 per hour have been achieved, as shown on Fig. [3U 

867 A detailed discussion of the machine performances after the implementation 

sea of the new interaction scheme can be found elsewhere [42|43] . The gain pro- 

869 vided by the new IR gets higher with the current products and the difference 

870 with respect to collisions with the crab sextupoles off can reach 50%. Figs. [35] 

871 and [HS] provide comparisons of the upgraded DA$NE with respect to past 
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Fig. 33. Typical example of a two hours running period at the upgraded DA&NE 
accelerator; these data have been taken between 00:52 and 02:52 on April 28 th 2009. 
The time evolution of the electron (blue) and positron (red) beam currents are shown 
on the top plot; the fill and cost pattern used to operate the machine is clearly visible. 
The middle plot stripcharts, among other quantities, the luminosities measured by 
the Bhabha calorimeter. The 'T1FREE' (green curve) is the raw luminosity directly 
computed from the detector trigger rate while the 'T2FARM' quantity (red dots) is 
background-subtracted. As expected, the latter is much less sensitive to backgrounds 
than the former: 'T1FREE' exceeds significantly 'T2FARM' when either beam gets 
injected while the two curves agree almost perfectly when the beam currents are both 
very low. With currents around 2 Amp for electrons and 1 Amp for positrons, the 
peak luminosity is routinely between 3.5 and 4.0 10 32 cm~ 2 s . Finally, the bottom 
plot shows the integrated 'T1FREE' and 'T2FARM' luminosities. Each step in the 
green curve corresponds to an injection period; after two hours the disagreement can 
reach 30% which explains why correcting the calorimeter rate was mandatory to get 
an accurate luminosity measurement. 
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Fig. 34. Two hours running period at high injection regime on December 21 rst 
2008; the integrated luminosity exceeds 1 pb^ 1 per hour while its instantaneous value 
oscillates around 3 10 32 cm _2 s _1 . This running mode cannot be maintained over long 
periods as it is not compatible with the SIDDHARTA operations. 



872 experiments (KLOE and FINUDA). The improvement is striking for all quan- 

873 tities measured: luminosity, specific luminosity and integrated luminosity. The 

874 saturation at high currents seen by the Bhabha calorimeter is also present in 

875 other monitoring data which indicates that the limitation does stem from a 

876 beam size blow up and not from some background-induced PMT saturation 

877 preventing accurate luminosity measurements. 
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Fig. 35. Comparison of the upgraded DA^NE performances (green dots) w.r.t. those 
achieved during the KLOE (blue and red dots) and FINUDA (yellow) runs. The left 
(right) plot shows the daily average luminosity (daily average specific luminosity) 
versus the product of the two beam currents. The improvement shown by the data 
taken on March 13 th 2009 is at least a factor 3 over the full range plotted on the 
x-axis. 



878 8.2 Bunch-by -hunch luminosity 



879 The Bhabha calorimeter timing resolution is good enough to separate the 

880 contributions from individual bunches, which is useful for machine studies 
sal (bunch filling pattern, etc.). Histograms of trigger TDC counts for the two 

882 pairs of back-to-back modules are shown in Fig. [37] for different accelerator 

883 configurations. 

884 • Top left plots: by-5 pattern (20+20 bunches). 

885 • Top right plots: by-3 pattern (33+33 bunches), 
sae • Bottom left plots: by-2 pattern (50+50 bunches). 
887 • Bottom right plots: full fill (100+100 bunches). 
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Fig. 36. Evolution of the daily and integrated luminosities for the FINUDA and 
SIDDHARTA runs of the DA&NE accelerator. The improvement provided by the 
crab waist compensation scheme is clear, both for the absolute performances and 
the rate of improvement. 



ass The TDC data are corrected by the RF timing signal: using the first bunch 

sag crossing as reference, one only keeps the time elapsed since this event which 

890 allows one to separate the contributions of each individual bunch to the trigger. 

891 The bunch spacing is found consistent for the different patterns: 2.6 ns in 

892 average. One can also note that the luminosity per bunch is quite flat for all 

893 bunch patterns and does not really depend on the bunch position in the train, 

894 as expected for a well-behaving machine. Obviously all fills are not as good 

895 as those exemplified in this section and so having the capability of making 

896 such diagnostics quickly helps a lot investigating issues which may involve the 

897 beam bunch patterns. 
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Fig. 37. Histograms of trigger TDC counts (1 count corresponds to 1.04 ns ) f or the 
two pairs of back-to-back calorimeter modules (M1-M4 and M2-M3) and different 
bunch patterns. Top left plots: by-5 pattern (20+20 bunches); top right plots: by-3 
pattern (33+33 bunches); bottom left plots: by -2 pattern (50+50 bunches); bottom 
right plots: full fill (100+100 bunches). These trigger TDC data are corrected by 
the RF timing signal to give the time elapsed since the last crossing of the first e~ 
and e + bunches; this allows one to see the contributions from individual bunches 
separately. 
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898 8.3 Touschek background 



899 Two machine runs have been dedicated to the measurement of the gamma 

900 monitor rates in order to compare the actual numbers with those predicted by 

901 a full tracking simulation. The latter simulates the Touschek scattering parti- 

902 cles from their generation to their loss in the beam pipe; secondary particles 

903 are then simulated and propagated through the IR up to the gamma monitors. 

904 The first run (February 3 rd 2009) was taken with two beams out of collision 

905 while the accelerator was operated in single beam mode for the second run 

906 (March 12 th 2009), during which the vertical beam sizes (and hence the cou- 

907 pling values) were varied using a skew quadrupole. For this second experiment 
90s a ten-bunch pattern was used for a total beam current of about 100 mAmp 

909 and the optics were the same as for the high luminosity runs. In addition to 

910 the gamma monitor spectra, the beam lifetimes were also measured during 
9n the March run. This allows a more complete test of the simulation as its pre- 

912 dictions regarding both the background rates and the beam lifetimes can be 

913 simultaneously compared with the DA$NE data. 

914 8.3.1 Gamma monitor rate 

915 The black triangles in Fig. show the variation of the gamma monitor rate 

916 versus threshold (in MeV), as obtained in simulation. The experimental setup 

917 makes the spectrum end around 440 MeV which is below the beam energy 

918 (510 MeV). Blue (red) bullets show the actual rates measured by the gamma 

919 monitor sensitive to photons emitted by positrons (electrons) for a few thresh- 

920 olds. In order to allow a direct comparison with the simulation, these numbers 

921 have been corrected to take into account the coupling between the beams. 
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922 Although the agreement between simulation and measurements is not perfect, 

923 both sets of numbers appear to be in the right ballpark - actual rates also 

924 depend on various machine conditions which are not easy to reproduce nor to 

925 quantify. 



Touschek rate 




10 =r 



100 200 300 400 500 

Energy (MeV) 

Fig. 38. Comparison between simulated (black triangles) and measured (blue dots for 
the electron beam, red ones for positrons) rates in the gamma monitors. The rates 
extracted from data have been corrected to take into account the beams coupling in 
order to allow a direct comparison with the simulation. 



926 8.3.2 Beam lifetimes 

927 The resulting simulated Touschek lifetime is as short as 840 s for a 0.5% beam 

928 coupling. This result is in agreement with measurements, as shown in Fig. [39j 

929 This plot compares the measured normalized lifetime of both beams with the 

930 computed one as a function of the square root of the effective beam coupling 

931 K e ff. Black markers refer to simulation while blue and red dots correspond to 

932 the measured electron and positron lifetimes, respectively. The plotted lifetime 

933 is normalized to the total current 

( i y /3 

934 ^normalized ''"measured I _, „„ a I \ ' ) 

\ 100 mAmp / 
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935 according to the scaling law r oc aia x a y /I where <j\ oc J 1 / 3 is the current- 

936 dependent bunch length. 

937 Larger and smaller markers refer to larger and smaller couplings respectively. 

938 The beam coupling 

939 K = e y /e x = (a v /a x ) 2 (3 x /(3y (8) 



940 is evaluated at the SLM using the measured transverse beam sizes and the 

941 ratio /3 y //3 x = 2.25, as indicated by the MAD optical model [16] . The effective 

942 vertical beam size used for the evaluation of the effective coupling K e ff takes 

943 into account its measurement resolution: 80 frni [40J. 
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Fig. 39. Comparison between predicted (black triangles) and measured (bullets, blue 
for the electron beam and red for the positron beam) beam lifetimes versus the ef- 
fective beam coupling K e ff - see text for details. The agreement between simulation 
and data is quite good. 
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944 9 Conclusions 



945 Various diagnostic detectors have been installed around the upgraded DA$NE 

946 IR to monitor and quantify the improvement of performances brought by the 

947 new crab waist collision scheme. After the design and building phases, the 

948 data taking started in February 2008 and is still going on in spite of a few 

949 shutdowns, in particular during the Summer 2008 break. 

950 The background-subtracted luminosity measured by the Bhabha calorimeter 

951 shows a significant improvement with respect to the previous IP design, valid 

952 for all currents and still increasing as the machine gets better tuned. This 

953 absolute measurement, based on a accurate GEANT simulation of the IR, 

954 has a systematics of about 13% and a negligible statistical error. In addition, 

955 gamma monitors allow online measurements of the background which has been 

956 found compatible with the predictions of simulation of the Touschek effect. 

957 Finally, the qualitative agreement between simulations and measurements is 

958 enforcing the validity of the studies currently ongoing to design the SuperB 

959 machine pQ. 
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